Fluctuations due to the nonlocal character of collisions 
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It is shown that the collision integral describing the nonlocal character of collisions leads to 
the same mean-field fluctuations as proposed by Boltzmann-Langevin pictures. It is argued that 
this appropriate collision integral contains the fluctuation-dissipation theorems in equilibrium itself 
and therefore there is no need to assume additionally stochasticity. This leads to tremendous 
simplifications in numerical simulation schemes. 
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The question how to describe sufficiently the fluctua- 
tions in dynamical systems of many interacting particles 
is as old as the discovery of the Brownian motion. Dif- 
ferent schemes have been developed in different branches 
of physics. Mostly additional stochasticity is assumed to 
account for such fluctuations. The goal of this paper is 
to show how the correct fluctuations can be described by 
a realistic collision scenario within the deterministic ki- 
netic theory without ad-hoc assumptions about stochas- 
ticity. Here we will briefly outline the attempts in nuclear 
physics having its counterparts in other fields, of course. 

It has been noticed in the beginning of the 80th that 
the time-dependent mean-field description (TDHF) of 
nuclear collisions cannot describe the experimental fluc- 
tuations of observables Q like mass, charge and momenta 
of emitted particles if the naive expectation 



wave functions 



AQ 2 (i) = Tr[Q 2 p n {t)} - [TrQp n {t)f 
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with the TDHF density matrix p n is applied. This 
deficiency has been cured by a variational approach 
|H EE El which leads to the propagation of fluc- 
tuations from to 

AQ 2 (h) = lim -Ltt^o.O) -p n (t ,v)] 2 (2) 

where p n (t,rj) evolves with the TDHF equa- 
tion but with the boundary condition p n (t, r/) — 
cxp (— ir)Q)p n (t) exp (irjQ). In this way the fluctuations 
are obtained by propagating back in time but they 
become explicitly dependent on the observables Q. It 
has been shown that the expression J5J leads besides the 
TDHF fluctuation (JIJ to an additional part which can be 
described by higher order diagrams in the interaction Q . 
The application of this procedure leads to a significant 
enhancement of the fluctuations Q, 0, Q ■ The fact that 
higher order diagrams are necessary to describe more 
appropriate fluctuations shows that the collisions are 
not described appropriately in usual Boltzmann (BUU) 
simulations. 

Alternatively there has been developed the time de- 
pendent generator coordinator method (TDGCM) 0,0] 
which expands the wave function $ in a set of TDHF 
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where the coefficients /jv are determined from minimizing 
the action. The TDGCM wave function leads to optimal 
fluctuations if the TDHF basic set 4>n(0) is chosen such 
that the variable Q becomes a generator of the path at 
finite time t. This TDGCM schema is equivalent to the 
above described Balian-Veneroni variational approach Q 
in the random phase approximation (RPA) limit. The 
advantage of the TDGCM schema is that it provides for 
optimal paths. However, both schemes are too limited 
for practical applications since one can handle only small 
sets of collective correlation channels. 

This practical limitation has led to the development 
of stochastic TDHF which approximates the time 
evolution of the N-particle density operator p at a small 
time steps ti by the diagonal elements of the expansion 
in TDHF density operators p n 

P (U) = £ W nPn (U) + W n P n >nU[p n >(U) ~ An(*i)]-(4) 
n nn' 

The transition probability is given by the matrix element 
of TDHF Slater states at time U 



= ^\{n\V\n')\ 2 S(E n -E' n ) 
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which leads to a Monte Carlo method of evolving an ini- 
tial state into an ensemble of Slater states with the prob- 
ability W m = 1-1] Pn'nU for m = n and W m = P mn 

n' 

for m ^ ti. 

The usually used Boltzmann collision terms with the 
inclusion of Pauli-blocking (BUU) cannot account for 
these fluctuations since the collisions are treated as ideal, 
i.e. local in space and time. Therefore there has been 
proposed another method of including more realistic fluc- 
tuations in the Boltzmann equation for the one-particle 
distribution / by adding a stochastic term 51 to the 
collision integral, called Boltzmann - Langevin picture 



collision integral, called r 

mm mfi Qua 13 



dt 



(1 - /)Win - fW mt + SI 
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where schematically the scattering - out and -in proba- 
bility of a phase space cell is W out /i n . This Boltzmann- 
Langevin equation can be formally derived from the 
stochastic TDHF equation if the one-particle reduced 
density is traced out. This equation has been applied 
for simulation of heavy-ion collisions quite frequently 

In these treatments the 
Langevin term in J^J is mostly assumed to be Markovian 



(81) = 0, (6I(t)6I(t')) = 2D6{t - t') 
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which leads with © to the equation of motion for the 
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variance er 
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with 1/t = Wout+Win- From the stochastic TDHF equa- 
tion (@J the form @ follows precisely neglecting fluctu- 
ations in the potential i.e. higher order diagrams |l0| . 
Again this is a hint that the collisions have to be treated 
more appropriately. Before suggesting a way to describe 
fluctuations more realistically, let us discuss some princi- 
ple problems of Boltzmann-Langevin approaches. 

We have in principle no reason to see the time evolu- 
tion of the density operator stochastically since the basic 
van-Neumann equation is deterministic and subsequent 
derived equations should be so as well. The ad-hoc as- 
sumption about stochasticity can mimic the numerical 
noise unavoidable in solving such equations and to a cer- 
tain extent higher-order correlations. Theoretically it is 
a problem since the collision integral emerges itself from 
averaging about small scale fluctuations [26j. Therefore 
it is ambiguous to divide correlations into an averaged 
collision integral and a stochastic term miming higher- 
order correlations. This has been sometimes motivated 
by the need of the fluctuation-dissipation theorem asso- 
ciating the collision integral with dissipation. In con- 
trast one should observe that the fluctuation-dissipation 
theorem emerges itself from appropriate collision inte- 
grals alone since they vanish in equilibrium. This can 
be seen best from ring summation approximation (RPA, 
GW, bubble^ leading to the Lenard-Balescu collision 
integral [27I Esf . The latter one vanishes exactly if the 
quantum fluctuation-dissipation theorem is fulfilled. In 
other words if the collision integral is derived appropri- 
ately it leads to an equilibrium with correct fluctuation- 
dissipation theorem and there is no need for additional 
stochastic terms. 

Introducing fluctuations the determination of / and 
the stochastic process 51 remains phenomenological in 
the sense that they account partially for such an appro- 
priate collision integral. This is a practical need if the 
appropriate collision integral is not solvable and one is 
restricted to Born (Boltzmann) approximation. Here in 
this letter we will show that the same mean fluctuations 
are generated from the deterministic but more realistic 



nonlocal extension of the Boltzmann collision integral 
[29l I30I l3ll l32]| . The advantage is that besides the micro- 
scopic foundation the latter one has already established 
a practical and fast numerical method. 

Let us return to the Boltzmann-Langevin equation. 
The fluctuation term derived from stochastic TDHF 
leads to the educated guess 



2D = (1 - (f))W in + (f)W out 
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which has suggested an interesting procedure to include 
fluctuations dynamically in BUU codes. From JHJ) one 
sees that the choice © reproduce the equilibrium vari- 
ance a 2 = /o(l — /o) an d that the deviation between 
the actual variance and the statistical value 8a 2 = a 2 — 
(/)(!-(/)) obeys 



dSa 2 2 c 9 
— — = — 6a . 
dt t 
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Since 5a 2 = is a solution of (|10() , the averaged value (/) 
can be parametrized as local equilibrium value such that 
at a given phase-space cell and time a small propagation 
will not deviate the variance from this result |2£|. Nu- 
merically the authors in [2(| used instead of the variance 
in distribution the density projection 
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which has been realized in each phase-space point of col- 
lisions. 

Besides the practical success of such descriptions in nu- 
merical solutions it remains the more principle question 
of the validity of the ad-hoc stochastic assumption. As 
we have pointed out this assumption partially cures the 
ideal collision scenario of space and time point-like parti- 
cles. We will follow now the other point of view and claim 
that if the collisions are described more realistically by 
nonlocal events then the fluctuations should be correctly 
induced by itself. 

The nonlocal extension of the Boltzmann equation has 
been given in |29ll30| and the finite duration and dynami- 
cal size of nucleons are calculated in |3lj | . The implemen- 
tation in BUU codes has allowed to describe experimen- 
tally noticeable effects j^U, I^J 0] • The prediction of a 
change in the reaction mechanisms j3^ has been nicely 
confirmed by proper scaling of experimental data [36j. 
The correlations by the nonlocal character of collisions 
are capable to describe long range order as typical for 
phase transitions [sflj . 

These nonlocal extensions are more realistic than the 
Boltzmann equation since they lead to the inclusion of 
two-particle correlations while in the Boltzmann equation 
these correlations are absent. We will show now that this 
more realistic description of the individual collision event 
induces also more realistic fluctuations which should be 
compared to the above described treatments. Especially 
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we will find that the nonlocal corrections lead to induced 
mean- field fluctuations of the form . 

The nonlocal quasiparticle kinetic equation for the dis- 
tribution function of particle a, f± = f a (k,r,t), derived 
within the non-equilibrium Green's function technique 
[2^, l30| has the form of a Boltzmann equation with the 
quasiparticle energy e± — e(k, r, t) 



dfi dei dfx 
dt dk dr 



ggi dfi 
dr dk 



•E 



dpdq 
(2tt) 5 ^ 7 
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hft{^-ff)^-ft) 
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and the spin-isospin etc degeneracy s. The superscripts 
± denote the signs of non-local corrections: f^ = 
f b {p,r±A 2 ,t), ff = f a (k-q±A K ,r±A 3 ,t±A t ), and 
f± = fb(p+q±AK,r±A4,t±A t ). For the out-scattering 
part of l|12l) both signs can be given equivalently [37J • The 
scattering measure is given by the modulus of the scatter- 
ing T-matrix V± = \T£\ 2 8{ei + e 2 ~ e 3 - e 4 ± 2A E ). All 
corrections, the A's, describing the non-local and non- 
instant collision are given by derivatives of the scattering 
phase shift <fi = Im \nT R (fl, k,p, q, t, r) 
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A 4 = - Nr + ^ 
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~dk 



dq 



E1+E2 
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The nonlocal kinetic equation i|12|) covers all quantum 
virial corrections on the binary level and conserves den- 
sity, momentum and energy including the corresponding 
two-particle correlated parts |30| . It requires no more 
computational power than solving the Boltzmann equa- 
tion mm. 

We will now derive the fluctuation contribution to the 
mean-field from this collision integral. Therefore we use 
the fact that the mean-field can be considered as zero- 
angle collisions, a collision which does not change mo- 
menta of particles but changes its phase. Summing up 
all these changes in phase gives just the mean-field po- 
tential due to the surrounding media. When deriving col- 
lision integrals one tries to share correlations in such a 
way that all mean-field like contributions are collected in 
the quasiparticle energies on the left side and all true col- 
lisions with finite transferred momenta are on the right 
side. In addition we will observe now that due to the 
nonlocal character of the collision there is a finite zero 
transfer momenta channel hidden in the nonlocal colli- 
sion integral which can be rewritten into the drift side of 
the kinetic equation and which gives exactly the fluctua- 
tions ijTHl . 



To show this we rewrite the energy-conserving 8- 
function in (|12fl as 
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where the channel q =/= represents the usually collision 
integral of l|12|) . The q — channel leads now to an 
additional part absent in usual local kinetic equations like 
the Boltzmann equation. To convince the reader about 
this novel observation let us rewrite the Pauli-blocking 
factors of 112[l for the q = channel according to (|14f> 



(1 - A) (I-/2") fsf* - hft (!- ft) (1 ~ ft) 

= / 2 ~(i-/ 2 ~)[/r-/i 



q=0 

(15) 



We see that in the case of local kinetic equations without 
delays this specific channel disappears since /{~ = f%. 
Therefore in contrast to the usual local kinetic equa- 
tions the nonlocal equation possesses a finite zero-angle 
channel in the collision integral which is of mean-field 
type since no energy or momenta is exchanged. The fact 
that retardation leads to an additional correction to the 
Bogoliubov-Hartree-Fock mean-field has been observed 
first within linear response in |3S| . 
From we obtain now 



dfi | ggi dfx 
dt dk dr 



dfi dfx 
dr dk 



(fr-fi) 



(27r) 5 fi 7 

(i-/ x ) (i-/ 2 -)/ 3 "/r - a^Ci-^) as) 



+s 



dpdq 



where 
1 



<^ASi^Gt - ^ M1 ~ (17) 



h J {2irhf 



m b 



Here we use the approximation of thermal averaged de- 
lays, (...}, which allows us to pull out the f^ term which 
contains the p-dependent shift under the integration of 
the 1/f term. This serves here for legibility and can 
be rendered exactly if one keeps the corresponding shift 
terms under the (p, q) integral of ifTTjl . 

The resulting equation (|16|l is a delay-differential equa- 
tion and has lead already to an interesting interplay be- 
tween stochastic bifurcations and relaxation due to iner- 
tia H3- We see from JTBJl and JTJJ that terms 1/f cx a 2 
of i|ll[) appears. In order to see that this additional term 
represents mean-field fluctuations let us expand 

(/i - a - )H a.!+a s §>4>^ + a4 (18 , 
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where in the last step we have replace the time derivative 
of / by the free drift motion « _<|lM + 
leading to the on-shell shifts 

A 3 = (A 3 



A* = (Ak + A^} = (^). (19) 

Now we can shift (11811 from the right to the left (drift) 
side of l |lti|l to obtain finally 




X [(1 -ft) (1 - / 2 -) / 3 -/ 4 - - ft ft (1- / 3 ± ) (1 - • (20) 

We see that an explicit fluctuating term ~ 1/f ~ cr^ (r, i) 
emerges to the mean field which has the form of l[lljl. 

Therefore we conclude that the variance of fluctuations 
from the Langevin-Boltzmann equation can be repro- 
duced in a deterministic way from the nonlocal extension 
of the Boltzmann equation. The stochastic treatment 
in numerical realizations can be considered therefore as 
a numerical trick to reproduce the correct determinis- 
tic fluctuations. The nonlocal collision scheme instead 
provides a first principle theory of such fluctuations. It 
gives probably the same practical results but with much 
less numerical effort since the collision scenario of the 
usual BUU code is modified only slight ly w ith no ad- 
ditional computational time required |32L l34j . Therefore 
we have demonstrated that the ad-hoc assumption about 
Langevin sources to the Boltzmann equation is unneces- 
sary if the collisions are treated non-locally. 

The clarifying discussions with P.G. Reinhard are 
gratefully acknowledged. 
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